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Density functional calculations are used to investigate the electronic structure of two-dimensional 5(i tantalum 
carbides with honeycomb-like lattice structures. We focus on changes in the low-energy bands near the Fermi 
level in pristine sheets of Ta^+iCn with increasing dimensionahty (n). We find that the Ta bd states dominate, 
but due to the extended nature of the wavef unctions, they are weakly correlated. The carbide sheets are 
prone to long range magnetic order, and we evaluate the stability of these states to enhanced electron-electron 
interactions, which we include using an onsite Hubbard U correction. Lastly, we find spin orbit interactions 
strongly renormalize the band structure for n = 2, but play a minor role in n = 1 and 3 sheets. 



Two-dimensional (2D) free-standing crystals exhibit 
a range of functional properties, mainly derived from 
the topology of their underlying lattice and enhanced 
electronic and magnetic effects due to reduced dimen- 
sionality. Spin-polarized edge states,^ for example, have 
been predicted for the well-studied 2D carbon material 
graphene,^ owing to the topological origin of its transport 
properties. A large external magnetic field, however, is 
required to realize the quantum Hall effect, and its spin 
degeneracy makes it difficult to manipulate. To overcome 
these challenges, experimental approaches have been de- 
veloped to induce magnetism by introducing transition 
metal adatoms^ and point defects^ on the surfaces. 

2D binary metal oxides and dichalcogenides, e.g. ZnO, 
BN, M0S2 also find widespread interest. New func- 
tionalities originate from the presence of more diverse 
chemistries.^'^ However, in most existing pristine 2D free 
standing materials, magnetic ordering is absent and the 
tunability of the electronic structure is limited to electro- 
static doping. A more promising avenue includes directly 
incorporating transition metals with multiple orbital de- 
grees of freedom and highly-correlated electrons into to 
the lattice. Magnetism, for example, was recently pre- 
dicted for VX2 (X=S, Se) monolayers.^ Alternatively, 
heavier bd transition metals with strong spin orbit cou- 
pling can be either deposited on the surface or directly 
integrated into the lattice of the 2D materials. Recent first 
principles calculations show that graphene decorated with 
bd transition metals can exhibit remarkable magnetic and 
topological transport properties.^ The magnetic coupling 
induced by such treatments on nonmagnetic 2D materials, 
however, is difficult to control in actual applications due 
to unintentional impurities and defects. 

In this Letter, we focus on novel low-dimensional hexag- 
onal materials, consisting of alternating layers of carbon 
and tantalum. These Ta-containing transition metal car- 
bides are part of a new group of 2D materials called 
"MXenes," recently synthesized by chemical exfoliation.^ 
Similar to the previously studied Ta-decorated graphene 
structures,^ these materials contain sheets of carbon in 
the inner layers and Ta atoms on the surface. In this case, 
however, the Ta layers are ordered, rendering them less 
susceptible to defects and more favorable for deliberate 



surface functionalization. Furthermore, these materials 
are derived from MAX phases, which are a large family 
of layered carbides and nitrides with the general formula 
Mn+iAXn^ where n = 1 • • • 3, M is an early transition 
metal, A is an A-group element (mostly from groups 13 
and 14), and X is carbon and/or nitrogen. Since the 
stacking sequence of the XMe octahedra in the hexagonal 
MAX phases depends on the stoichiometry, the MXene 
sheets have the advantage that dimensionality controls 
both the system size and the symmetry between the two 
surfaces (Fig. la). 

The 2D MXenes have the general formula M^+iX^, and 
they crystallize in sheets containing 1, 2, or 3 layers of 
XMe octahedra depending on n. Their recent synthesis^ 
has led to a number of theoretical studies, ^^^^ mainly 
focusing on the Ti-containing phases. However, MXenes 
in the Ta-C system have not yet been well studied. 

Motivated by the structural fiexibility, possible en- 
hanced electron-electron interactions, and the strong 
spin-orbit coupling parameter for Ta (Cd = 1970 cm"-*^),-^^ 
we use ah initio electronic structure calculations to in- 
vestigate the effect of dimensionality (n) and electron 
correlations on the band structure and magnetic ordering 
in Tan+iCn, n = 1 • • • 3. 




FIG. 1. The two-dimensional Tan+iCn (n = 1,2,3) sheets (a) 
all possess hexagonal symmetry. The lattice constants, a, and 
select interatomic Ta-C distances are highlighted (cf. Table I), 
(b) LDA electronic band structures for each compound along 
the path r(0,0,0) ^ M(|,0,0) ^ K[\, |,0). 
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TABLE I. Summary of lattice parameters, a, and Ta-C bond 
lengths, d, in A with respect to dimensionality, n, obtained 
with LDA and PBEsol functionals. Embolden values corre- 
spond to experimental data taken from Ref. 19. 
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First-principles density functional calculations are per- 
formed using the Vienna Ab initio Simulation Package 
(VASP)/^ with a plane wave cutoff of 500 eV and the 
projector- augmented wave method (PAW)^^ to treat the 
interaction between the core and valence electrons; we 
treat the Ta 5p electrons as valence electrons. Reciprocal 
space integrations are performed using a 9 x 9 x 2 /c-point 
mesh. We investigate the effects of electron-electron 
interactions by using both the local (spin) density ap- 
proximation [L(S)DA)] and the improved generalized gra- 
dient approximation (GGA) of Perdew-Burke-Eruzerhof 
(PBEsol) for solids with the rotationally invariant Hub- 
bard U correction (+/7) of Liechtenstein et al.^^^ keeping 
J = 0.5 eV throughout. The high measured conductivity 
of Ta4C3,^^ and the robustly metalhc electronic structures 
of bulk TaC^^ and the MAX phases, suggest that any 
electron correlation effects should be weak. We, therefore, 
anticipate the LDA and PBEsol functionals to provide 
an adequate description of the electronic and magnetic 
properties of these materials. 

The Ta^+iCn unit cells used in our calculations con- 
tain two symmetry equivalent free-standing sheets that 
are separated by 11-13 A of vacuum. We obtained the 
equilibrium structures at the LDA and PBEsol level by 
minimization of the total energy computed for a range 
of a lattice parameters, performing a full relaxation of 
the atomic positions along the c-direction until the forces 
were converged below a tolerance of 5 meV A~^. 

Table I contains the ground state atomic structure de- 
scriptors obtained with the LDA and PBEsol functionals. 
The LDA functional, for all values of n, predicts equi- 
librium lattice parameters that are smaller than those 
obtained with PBEsol; nonetheless, both functionals are 
in good agreement, within < 1% of each other and the 
available experimental data^^ for Ta4C3. We also sum- 
marize the interatomic distances between the different 
Ta and C atoms corresponding to the sites labeled in 
Fig. 1(a). For all 2D sheets explored, the Ta atoms at the 
surface layer have shorter Ta-C bonds than those in the 
center of the sheet, i.e. (i(Tai-C) < (i(Tan-C). 

While there are slight differences in the lattice dimen- 
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FIG. 2. Site-decomposed partial DOS computed within the 
LDA for the surface (surf) and inner (in) layer Ta sites in (a) 
Ta2C, (b) Ta3C2, and (c) Ta4C3. Note that the contribution 
from the C atoms to the DOS within this energy window 
is small and therefore not shown. The spatial distribution 
of the partial charge density from 0.2 eV up to the Fermi 
level (shaded region) is also shown in (d), (e), and (f) for 
Ta2C, Ta3C2, and Ta4C3, respectively. The 2D contours are 
projections on the (1120) plane. 

sions and atomic positions, the average interatomic dis- 
tances and lattice parameters are similar between the 
three stoichiometrics. This suggests that any differences 
in electronic structure should originate from either the 
stacking of the octahedra or the number of occupied Ta 
d bands available, which depends on the ratio of Ta to C 
atoms. To explore these possible differences, we begin by 
computing the electronic band structures with the LDA 
functional. For n = 1, we find two dispersive Ta d-bands 
crossing the Fermi level (Ep) along F — M [Fig. 1(b), 
left] . Other than the points near F where these two bands 
cross Ef^ the Fermi level falls in the band gap present 
along most of the sampled path in the Brillouin zone (BZ). 
These two free-electron pockets centered at F are similar 
to those found in 3D metals despite the 2D nature of the 
MXene sheet. In contrast, we find in Ta3C2 (n = 2) a 
single Ta band that crosses the Fermi level along F — M 
with nearly linear dispersion [Fig. 1(b), center]. Due to 
the symmetry of the BZ, there is another linear band 
crossing of this type along K— F. At the M-point, the 
band just below Ejr, corresponding to a Ta d orbital, is 
fully occupied. The band structure for Ta4C3 (n = 3), 
on the other hand, is strongly metallic; three Ta d bands 
cross Ef [Fig. 1(b), right]. The lower two bands rise 
above Ef around the M-point and again cross the Fermi 
level at the K-point. Indeed the site-decomposed partial 
densities-of-states (DOS) confirm that the region near Ef 
is largely controlled by the Ta d states [Figs. 2(a-c)]. 

Figs. 2(d-f) show the spatial distribution of the elec- 
trons within 0.2 eV of ^i?. We find a strong dependence 
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FIG. 3. The L(S)DA+[/ energy differences, A^, between the 
magnetic and non-spin polarized (NSP) states with U for (a) 
Ta2C, (b) Ta3C2, and (c) Ta4C2. Energies given per Ta atom. 
The points at [/ = eV are obtained without the plus Hubbard 
U method. Schematics (inset) illustrate the spin ordering of 
the Ta atoms with ferromagnetic (FM), ferrimagnetic (FiM), 
or antiferromagnetic (AFM) order; the arrows represent the 
relative magnitudes and directions of spin, and the open circles 
indicate a Ta atom with no magnetic moment. 

on dimensionality for the charge distribution. Intrigu- 
ingly, this spatial distribution about the Ta site in n = 1, 
viz Ta2C [Fig. 2(d)], and the inner Ta atom in n = 3, 
viz Ta4C3 [Fig. 2(f)], share similar features — the charge 
around the atom is distributed into six lobes, three above 
and three below the Ta atom. The inner Ta atom in n = 2 
(Ta3C2) [Fig. 2(d)] shows strikingly different behavior, 
with a small distribution of charge collected above and 
below the Ta atom, aligned parallel to the c-axis. 

To understand the atomic-scale origin of these features, 
we examine more closely the crystal structures of each 
sheet. Ta2C (n = 1) consists of a single Ta-C octahedron 
with stacking sequence AcB, where uppercase letters de- 
note the Ta atom and lower case letters represent the 
stacking of carbon atoms. Ta4C3 (n = 3), therefore, is 
obtained as a Ta2C layer, AcB, with an extra Ta atom on 
each surface - that is, Cb[AcB]aC. Note that this is the 
only structure in which two Ta atoms on either surface 
have the same stacking. On the other hand, Ta3C2 {n = 2) 
has a stacking sequence of AcBaC, so the symmetry of the 
AcB layer is broken due to the odd number of layers. This 
suggests that the differences in the charge distribution we 
find in Figs. 2(d-f) are largely governed by the stacking 
sequence of the CTae units — a degree of freedom unique 
to the MXene phases. Such differences in the electron 
distribution near Ep due to the asymmetric stacking are 
also visible in the shape of the DOS [shaded region in Fig. 
2(b)]. The charge is concentrated at the surface for Ta3C2 



TABLE 11. Summary of the stable spin polarized ground 
states for the Tan+iCn MXene phases using different exchange- 
correlation functionals with and without a Hubbard U correc- 
tion. Notations (schematics) for the magnetic states are given 
in the caption (insets) to Fig. 3. 
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and Ta4C3, consistent with the higher partial density of 
states for the Ta surface atoms [Figs. 2(b-c)]. 

The large number of Ta d-states at Ef and the sensi- 
tivity of the electronic structure to the sheet dimensional- 
ity suggest the possibility of stable long range magnetic 
spin configurations. We therefore performed a series of 
spin-polarized calculations with two different starting con- 
figurations, corresponding to ferromagnetic (FM) and 
antiferromagnetic (AFM) spin order on the Ta sites, to 
systematically explore the possible magnetic orders with 
respect to n. We carried out unconstrained-spin density 
calculations to find the stable magnetic ordering (Fig. 3) 
and compared the total energy of those states to that 
of the non-spin-polarized case. Given the limited abil- 
ity of DFT to fully capture correlation effects, including 
transition metal ions with partially filled d shells, we now 
add the Hubbard U correction to the standard PBEsol 
(PBEsol+/7) and LSDA (LSDA+/7) functionals. Since 
it was recently suggested that Ta-containing MAX phases 
are weakly correlated, we explored U < 3.0 eV. 

Figure 3 shows the change in total energy computed 
with the LSDA+/7 functional compared to the non-spin- 
polarized (NSP) case as a function of U. In all cases, the 
magnetically ordered configurations are lower in energy 
than the NSP state. As the dimensionality increases, the 
AFM case becomes more stable, likely due to the fact 
that the Ta atoms on each surface are more strongly spin- 
polarized. Our main results are summarized in Table II, 
where the values in parentheses specify the U values above 
which the specific magnetic ordering becomes stable. For 
Ta2C and Ta3C2, both LSDA+/7 and PBEsol+/7 calcula- 
tions predict a FM configuration, whereas an AFM order- 
ing is predicted for Ta4C3, but only with the LSDA-ft/ 
exchange- correlation functional (Table II). Note that in 
Ta4C3 no magnetic configurations were found to be stable 
using either PBEsol and PBEsol+[/. 

In general, a larger Hubbard U value leads to a more 
stable ferromagnetic (or ferrimagnetic) configuration, es- 
pecially for Ta3C2 [Fig. 3(b)]. In all cases, the surface Ta 
atoms are spin polarized in the same direction for FM 
and FiM order and in opposite directions for AFM. The 
Ta atoms in the inner layers are weakly spin-polarized in 
the direction opposite to the surface Ta atoms in both 
the FiM and AFM configurations, with the exception of 
the AFM configuration of Ta3C2, which is constrained by 
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FIG. 4. LDA band structure for Ta3C2 without (black, broken 
lines) and with (red, solid lines) spin-orbit coupling (SOC). 

symmetry. The addition of U generally has a small effect 
on the electronic structure, leading to a slight shift of the 
Ta d-bands uniformly to higher energies. In Ta3C2, whose 
band structure is controlled largely by a single band that 
is partially filled, the addition of correlation effects drives 
a transition from a state with a single band crossing Ef 
to a multi-band metal with a complex Fermi surface. 

Having established that the long-range magnetic order 
due to the spin exchange splitting of the Ta d states at 
Ef strongly alters the orbital filling, we now evaluate the 
effect of spin-orbit coupling (SOC) on the band structure 
with respect to n. Here we find that SOC splits the band, 
manifesting as a shift of the Ta d-bands to higher energy 
levels in all cases; it has the strongest effect on Ta3C2 
(Fig. 4), where the band structure changes from a nearly 
single band at the Fermi level to one with multiple band 
crossings throughout the Brillouin zone. The addition of 
SOC destroys the single band crossing in manner similar 
to that found for a finite U value. The spin-orbit split- 
ting at the top of the valence band (A5'o=148 meV) is 
approximately half of that observed in GaAs (A5'o=342 
meV),^'^ but more than 3 orders of magnitude larger than 
graphene {A so =~ 0.05 meV).^^ The result is that a sin- 
gle band n-type metal with linear dispersion transforms 
into a multi-band p-type conductor with hole pockets, 
which should exhibit a strongly temperature dependent 
thermopower. 

In summary, we have shown that the Ta-based bd elec- 
tronic structure is highly sensitive to dimensionality. All 
explored phases exhibit correlation stabilized magnetic 
order that is not found in the bulk MAX phase struc- 
tures. The LSDA+/7 method stabilizes the ferromagnetic 
ordering in the case of n = 1,2, and for n = 2 spin- 
orbit couphng renormalizes the electronic band structure 
with a transition from a nearly filled single band to a 
metal. In these 2D MXenes, the electronic structure is 
controlled by the stacking of the XTae octahedra and the 
states derived from the surface Ta atoms. Tailoring the 
electronic structure could therefore be achieved through 
end group functionalization of the surfaces of the MXene 
sheets. This opens up possibilities for engineering a new 
class of tunable functional 2D materials. We conjecture 
that one could maintain a graphene-like band structure 



like that observed in Ta3C2 by through epitaxial strain 
engineering. Overall, the Ta-containing graphene-like car- 
bides show great promise as functional 2D materials that 
can be synthesized in different dimensionalities, leading 
to a range of stacking sequences and stoichiometrics that 
offer a variety of electronic and magnetic behaviors. 
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